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ABSTRACT 

Using the data taken from Tibet II High Density (HD) Array (1997 February- 
1999 September) and Tibet-Ill array (1999 November-2005 November), our pre- 
vious northern sky survey for TeV 7— ray point sources has now been updated by 
a factor of 2.8 improved statistics. From 0.0° to 60.0° in declination (Dec) range, 
no new TeV 7— ray point sources with sufficiently high significance were identi- 
fied while the well-known Crab Nebula and Mrk421 remain to be the brightest 
TeV 7— ray sources within the field of view of the Tibet air shower array. Based 
on the currently available data and at the 90% confidence level (C.L.), the fiux 
upper limits for different power law index assumption are re-derived, which are 
approximately improved by 1.7 times as compared with our previous reported 
limits. 
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1. Introduction 

The development of TeV 7— ray observations has experienced a revolutionary progress'^'^l 
since the finish of our northern sky survey work'^l (hereafter Paper I). For example. High 
Energy Stereoscopic System (HESS) experiment alone has discovered more than 40 new 
7— ray sources in southern hemisphere with unprecedented angular resolution and sensi- 
tivity. Together with other sensitive Imaging Air Cherenkov Telescopes (lACTs), such as 
Major Atmospheric Gamma Imaging Cherenkov Telescope (MAGIC), Collaboration of Aus- 
tralia and Nippon (Japan) for a GAmma Ray Observatory in the Outback (CANGAROO), 
Very Energetic Radiation Imaging Telescope Array System (VERITAS), more than 50 new 
sources have been discovered in the past several years and the number as well as the di- 
versity of TeV 7— ray sources have been increasing. Spatial and temporal information of 
these sources are now available with high accuracy and this makes it possible for further 
studies on acceleration of high-energy cosmic rays, relativistic astrophysics, as well as quan- 
tum gravity theory and so forth. To demonstrate the advantage of its wide field of view 
and high duty cycle, the Tibet Air Shower Array experiment performed a high-precision 
measurement on the two-dimensional (2D) anisotropy of cosmic rays in the energy range of 
a few to several hundred TeV and discovered a fairly compact new anisotropic component 
in the direction of Cygnus region^ . Furthermore, MILAGRO experiment has discovered an 
extended 7— ray source in the direction of Cygnus region^^l and a few more other sources in 
the Galactic plane'^l, in addition to the diffuse 7— ray emission from the Galactic plane'^l. 
While some of the MILAGRO sources were confirmed or supported by the Tibet Air Shower 
experiment [^'^'^1, it would be extremely interesting and important for the Tibet Air Shower 
experiment to systematically update its northern sky survey with a much larger data sample 
currently available. 

2. Tibet Air Shower Array Experiment and Observations 

The Tibet air shower array experiment has been successfully carried out at Yangbajing 
Cosmic Ray Station (90.522°E, 30.102°N) in Tibet, China, since 1990, at an altitude of 
4300m above sea level. Having been upgraded several timest^^'^^'^^^, the Tibet HD and III 
arrays have identical structures except the array size and shape. A 0.5 cm thick lead plate 
was later placed on top of each counter to improve fast-timing (FT) data by converting 7 rays 
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into electron-positron pairs. The angular resolution was first estimated from Monte Carlo 
(MC) simulations and then confirmed experimentally by observing the Moon shadow to be 
about 0.9° in the energy range above 3 TeV. The data used in this analysis were collected 
by running the Tibet HD array for 555.9 live days from 1997 February to 1999 September 
and the Tibet III array for 1318.9 live days from 1999 November to 2005 November. The 
events are selected by imposing five criteria on the reconstructed data: (1) Each shower event 
should fire four or more FT detectors recording 1.25 or more particles. (2) The estimated 
shower center location should be inside the detector array. (3) ^ pft should be larger than 
15, where X] Pft is the sum of the number of particles per square meter detected in each 
detector. (4) The zenith angle of the incident direction should be smaller than 40°. (5) The 
residual error in direction reconstruction should be less than 1.0m. After applying these 
cuts and a data quality controU, about 2.0 x 10^° shower events were available for our data 
analysis here. 



Based on the successful analysis of Paper I and for simplicity, Method II (i.e., the all- 
distance "equi-zenith angle" method) was adopted to construct the 2D cosmic ray intensity 
map with pixels in the size of 0.1° x 0.1° in equatorial coordinate. The idea of this method 
is that at any moment, for all directions, if we scale down (or up) the number of observed 
events by dividing them by their relative cosmic ray intensity, then those scaled numbers 
of events in a zenith angle belt should be equal anywhere in the sense of statistics. A 
function can be built accordingly, the relative intensity of cosmic rays I {R. A., Dec) and its 
error AI(R.A., Dec) in each direction can be solved by minimizing the function. It is 
worth mentioning that source information in a 0.1° x 0.1° bin had been included in the 
intensity I [R. A., Dec). For details of this method, the reader is referred to Paper I. To 
remove the large-scale cosmic ray anisotropy and to keep the local event excess structure 
which is due to the 7— ray emission, we use the similar subtraction procedure as in Paper 
I when parameterizing the projected intensity distribution along the right ascension (R.A.) 
direction for any Dec belt. After subtracting the anisotropy, we can obtain the relative 
intensity of cosmic rays Icorr{R-A., Dec) and its error AIcorr{R-A., Dec). The number of 
excess events and their uncertainties in cell (R.A.,Dec) can be calculated as 



3. Data Analysis 



Ns{R.A.,Dec) = [Icorr{R-A., Dec) - l]Nobs{R-A.,Dec)/I, 



corr 



{R.A., Dec) 



(1) 



AN,{R.A.,Dec) = Ahorr{R-A., Dec)Nobs{R-A., Dec) / he 
where Noi,s{R-A., Dec) is the number of events in an on-source bin. 



corr 



{R.A.,Dec) 



(2) 
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Given the angular resolution of the Tibet air shower array, events are summed up from 

a cone with an axis pointing to the source direction, and the half-opening angle is set as 
0.9° (for E > 3TeV)[^3l or 0.4° (for E > lOTeV). All celestial cells with their centers located 
inside the cone contribute to the number of events as well as its uncertainty. Finally, the 
significance for an on-source window centered at the cell {R.A.on, DeCon) can be calculated 
by 

^ _ E(R.A.,Dec)econe{[^corr{R.A.,Dec) - l]No,,{R.A., Dec) / IcorrjR.A., Dec)} 

b\R.A.gri,i J^GCon) — , 

^JT.iR.A.,Dec)econe[^Icorr{R■A., Dec)Nobs{R.A., Dec) / Icorr{R.A., Dec)] 

(3) 

The systematic uncertainty for the significance value due to the above subtraction procedure 
on large-scale anisotropy is estimated to be 0.2a by adjusting the bin size and the smoothing 
parameters. 



4. Results and Conclusions 

Distribution of significance for all bins in the surveyed sky are shown in Fig.l. It 
agrees very well with a normal distribution on the negative side, indicating that systematic 
effects are well under control. The positive side contains more high-significance entries 
than those expected from pure statistical fluctuations, and they are related to two well- 
known TeV 7— ray sources, namely the Crab Nabula and Mrk421. After removing their 
contributions, in such a way that those cells within 2° regions around the Crab Nebula and 
Mrk421 are excluded, we get the dash-dotted histogram as shown in Fig.l, consistent with 
the expectations from random background fluctuations. 

Fig.l The significance map is shown here. The solid curve is derived from all cells defined 
in the analysis. The dash-dotted histogram excludes cells close to the Crab Nebula or 
Mrk421. The dashed line represents the best fit of a Gaussian curve to the data, its mean 
is -0.002 ± 0.01 and standard deviation is 1.013 ± 0.005. 

For our updated sky survey, the information of flve candidates for possible 7— ray sources, 
each with an excess of greater than 4.5(7, arc summarized in Table 1. Only the pixel with 
the highest signiflcance from each independent direction is listed. 

As can be seen in Table 1, the hst includes two estabhshed sources Crab Nebula and 
Mrk421 which remain to be the brightest TeV 7— ray sources in the northern sky. Compared 
with Paper I, the significance of the Crab Nebula is increased from 5.0a to 7.1a at the highest 
signiflcance position (from 4.1cr to 6.0a at the nominal position of the Crab Nebula, which 
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Signiflcance 



T able 1: Candidate locations of all directions with an excess greater than 4.5 (7 



No. 


R.A. 


Dec 


NoN 


Nqff 


Ns 


ANs 


^pretrials 


1 


57.95 


53.25 


2405072.8 


2397926.7 


7146.1 


1548.5 


4.6 


2 


70.55 


11.35 


2306840.6 


2299785.4 


7055.2 


1516.5 


4.7 


3" 


83.75 


21.95 


3078848.1 


3066434.9 


12413.3 


1751.1 


7.1 


4 


89.45 


30.05 


3359526.5 


3350799.7 


8726.8 


1830.5 


4.8 


5" 


166.25 


38.25 


3301780.3 


3292945.8 


8834.4 


1814.6 


4.9 



Table 1 The columns arc (from left to right) sequence of prominent direction, R.A.(.J2000), 
Dec(J2000), number of measured events in on-source window (Nqn), number of background 
events {Nqpf), event number excess in on-source window ( Ns=Non-]^off), uncertainty 
on the event number excess (ANg), and the significance Spretriais of deviation Nqn from 

NoFF- 

Notes: R.A. and Dec columns are due to the way we divide the bin in the analyses; 3° 

The Crab Nebula and 5^ Mrk421. 

is consistent with the expected enhancement from 3.7o" to 6.2a within error bars according 
to the changing statistics); while the significance of Mrk421 is dropped somewhat mainly 
due to the fact that Mrk421 is not a stable source; it happens to be in a high state with 
the data used in Paper I but remains less active'^^' in the succeeding period when the data 
are used in the current analysis. As for the other two candidates in Table 1, No. 2 and No.4 
have also appeared in Paper I; however their significance values are slightly decreased after 
we have included more data into this analysis. There are still two other sources located at 
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(70.45°, 18.05°) and (221.75°, 32.75°) which passed in Paper I but failed this time, and the 
current analysis selects a new candidate at position (57.95°, 53.25°) with a significance value 
just above 4.5(7. It should be noted that the above— discussed phenomena, except the Crab 
Nebula and Mrk421, are probably due to the background fluctuations and also possibly due 
to their intrinsic unstable features; the conclusive results will rely on a further data analysis. 
In summary, compared with Paper I, the number of hot spots reduces from 4 to 3 (only 
for the candidates determined from Method II and by excluding the two known sources: 
the Crab Nebula and Mrk421). Both are consistent with the expectation from statistic 
fluctuation. With 200 Toy MC experiments, the numbers of hot spots(each satisfles 4.5(7 
requirement) are obtained for each experiment, the probabilities to observe no less 4 and 3 are 
found to be 8% and 26% respectively. In addition, the locations of candidates are somewhat 
different between the two observations. As for the difference, it agrees with the pointing 
accuracy of the Tibet AS7 array. Taking Crab Nebula (0.4° position difference between two 
observations) as an example, we estimate the probability with MC experiments and find 40% 
of the MC experiments have a position difference no less than 0.4°. Nevertheless, given the 
large number of trials, the significance values from all directions other than the Crab Nebula 
and Mrk421 are not high enough to definitely claim any existence of a new point source, 
although they will become clearer with the future improved statistics of observational data 
or can be interesting regions for further follow-up observations with more sensitive lACTs. 

It is worth mentioning that the two MILAGRO newly reported TeV sources, MGRO 
J1908+06161 and MGRO J2019+37151 not fisted in Table 1 due to their lower significance, 
had been our two candidates with only marginal yet persistent event excess. In Paper I, we 
found 4.8c7 on (286.65°, 5.55°), 0.4° angular separation from MGRO J1908+06. However 

the significance value is 4.3(7 in the current analysis, not scaled up with the increase of 
statistics but consistent with the expectation based on the flux measured by MILAGRO'^' and 
HESS'^"''^'. Another interesting point source, close to an extended source MGRO J2019+37, 
has been discussed several times in Tibet AS7 papcrs''''^®'^'''. Our dedicated analysis has 
reported a preliminary 5.8(7 excess in (304°, 36.1°)[^1 and this result should be regarded as 
a confirmation of the MILAGRO's discoveryPl. While less sensitive to such an extended 
source, the current point source search analysis still finds a 4.0(7 excess in the direction of 
(303.25°, 35.95°). 

Based on the above analysis we know that the significance from all directions other 
than the Crab Nebula and Mrk421 are not high enough to definitely claim any existence of 
a new point source, we set a 90% C.L. upper fiux limit for all directions in the sky, except 
at the positions of the Crab Nebula and Mrk421. The prescription of Helene'^^' is used to 
calculate the upper limits of the number of signal events at the 90% C.L. for energies higher 
than 3TeV and lOTeV from each region of the northern sky. Then the effective detection 
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area of the Tibet air shower array is evaluated by full MC simulation assuming a Crab-like 
7-ray spectrum ^"2-6 for a set of Dec values (0.0°, 10.0°, 20.0°, 30.0°, 40.0°, 50.0°, and 
60.0°) and interpolated to other Dec values between 0.0° and 60.0°. Taking into account the 
live time, the newly derived 90% CL. average flux upper limit along the R.A. direction as a 
function of Dec is shown in Fig. 2(a), which is (0.8 ~ 1.9) x 10~^^ cm~^ s~^ for E > 3TeV 
and (1.3 ~ 2.5) x 10~^^ cm~^ s~^ for E > lOTeV respectively. Additionally, since the 
response of the Tibet air shower array is energy dependent, the flux upper limits obtained 
from these data are dependent on the energy spectra of the possible sources of TeV gamma 
rays. The same procedure is applied to the cases of other power-law indices for energy above 
3TeV and lOTeV, the corresponding average flux limits can be found in Fig.2(b), which are 
(0.8 ~ 2.2) X 10^12 s"^ for E > 3TeV and (1.2 ~ 3.0) x lO'^^ ^.^^-2 g-i ^ > lOTeV 

respectively and have approximately 1.7 improvement compared with the reported limits in 
Paper I. These limits are well consistent with the fact that the majority of 7— ray sources 
discovered in recent years have integrated fluxes less than 10% of those of the Crab Nebula 
at ITeV energyl^^l. 



90% CL. Upper Limit ( E>3 TeV , E>10 TeV ) 
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Fig. 2 R.A. direction-averaged 90% CL. upper limit on the integral flux above 3TeV 
and lOTeV. (a) for a Crab-like point source, i.e., with an energy spectrum of E^^-^; 
(b) for different indices of power-law spectra. 



In conclusion, we performed an updated northern sky survey for the TeV 7— ray point 
sources in a Dec band between 0.0° and 60.0° using about eight-year data obtained from 
February 1997 to November 2005 by the Tibet air shower array. The significance except Crab 
and Mrk421 is not high enough to definitely claim any existence of new sources. Accordingly, 
more stringent 90% CL. flux upper limits than the ones in Paper I are set from the rest of 
positions based on the assumption that candidate point sources have power-law spectra with 
indices varying from 2.0 to 3.0. In the near future, we will add a large muon detector array 
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under the Tibet air shower array for the purpose of increasing its 7— ray sensitivity in the 
100 TeV energy region (10-1000 TeV) by discriminating between 7 rays and the cosmic— ray 
hadrons'^^l According to a full MC simulation, flux sensitivity of this new project will be an 
order or more better than the present one in the lOOTeV regionl^^l Approximately 10 new 
sources are expected to be discovered and we will be able to measure the cutoff energies of 
known and unknown sources which are potential origins of Galactic cosmic rays. 

The collaborative experiment of the Tibet Air Shower Arrays has been performed under 
the auspices of the Ministry of Science and Technology of China and the Ministry of Foreign 
Affairs of Japan. This work is supported in part by Grants-in-Aid for Scientific Research on 
Priority Areas (712) (MEXT), by the Japan Society for the Promotion of Science (JSPS), 
by the National Natural Science Foundation of China (in part by NSFC Grants 10675134 
and 10533020) and by the Chinese Academy of Sciences. 
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